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v_ Dice Probability 


Working out the most probable round the game finishes on. 


This game consists of two teams. 


¢ Team A has one unit with 5 health. 


¢ Team B has infinite units, each unit has 1 health. 


The game follows a turn based system. 


Every turn there are two phases. 


¢ Phase 1: Duel Roll 
e Phase 2: Attack Roll 


Phase 1: Duel Roll 


¢ Team A rolls 4 dice 


¢ Team B rolls 8 dice 


Which ever team has the highest dice value wins the phase. On a draw, Team A win. 


e.g. 


Team A dice: 4, 2, 3, 5 
Team B Dice: 3, 5, 2, 4, 3, 1, 4, 3 


Team A wins becuase both teams rolled a 5, and Team A wins on draws. 


Phase 2: Attack Roll 


Which ever team wins Phase 1 gets to attack. 

If Team B win Phase 1 and get to attack, they make 16 rolls (2 * 8). 
To take health off they must roll two 6's in a row. 

Each success removes one health from the unit in Team A. 


Saving Roll: If the unit in Team A has one health, and an attack is succesful against it, then Team A may roll one dice, on a 2 or above the 
attack becomes unsuccesfull. 


Team A also has 4 skill points, this allows Team A to increase the score of this dice roll 4 times, this skill point will only ever be used to 
increase a Saving Roll if tyhe saving roll is a 1. 


e.g. 


Team A unit has 2 health left. 

Team B wins Phase 1. 

Team B makes 16 rolls. 

5 of these rolls are a 6. 

These 5 dice are re-rolled. 

2 of these re-rolled dice land on 6. 

Team A loses 1 health due to one of these dices. 

Team A loses another 1 health due to the other dice, however, Team A can make a saving roll. 
Team A rolls one dice. 

The result is a 1. 

Team A has 2 skill points left, so use one of these skill points to increase the result of 1 into a 2. 


Team A does not lose its last health. 
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End Turn 


¢ Now the Attack Roll has conducted, go back to the duel roll and start the process again. 
¢ Continue until Team A loses all health. 


Phase 2: Attack Roll (Team A wins Duel Roll) 


¢ Team A wins Duel Roll. 

¢ Team A roll 8 dice. 

e Every dice over 3 "kills" a unit on Team B. 

e Units on Team B are replaced after every Attack Phase infinitely. 


Goal of this excercise 


¢ Item 1: Probability of Team A or Team B Winning Phase 1. 

e |tem 2: Determine the probability of Team A wins Phase 1 in x amount of rounds out of j total rounds. 
e Item 3: Determine the probability of Taking 4 Health off Team A with no roll save. 

e Item 4: Determine the probability of Taking 4 Health off Team A with roll saves. 

e Item 5: Determine the probability of taking all Health off off Team A. 

e Item 6: Determine the probability of killing j amount of Team B units. 


v Item 1: Probability of Team A or Team B Winning Phase 1. 


¢ This is the number of combinations of dice results that can occur divided by the total number of combinations available. 


Imagine a game where instead of the above dice quantities, Team A can roll one dice, and Team B can roll 2 dice. The following outcomes are 
possible for Team A to win. 


Team A | Team B 

a | (1) and (1) 

2 | (1 or 2) and (1 or 2) 

3 | (1 or 2 or 3) and (1 or 2 or 3) 

4 | (1 or 2 or 3 or 4) and (1 or 2 or 3 or 4) 

5 | (1 or 2 or 3 or 4 or 5) and (1 or 2 or 3 or 4 or 5) 

6 | (1 or 2 or 3 or 4 or 5 or 6) and (1 or 2 or 3 or 4 or 5 or 6) 


However, the results of Team B actually can be any combinations of the results: 


Team A | Team B 

1 | (1 and 1) 

2 | (1 and 1) or (1 and 2) or (2 and 1) or (2 and 2) 
etc. 


Count the number of different combinations that can occur: 


Team A | Team B 

4 | 1 combination 
2 | 4 combinations 
3 | 9 combinations 
4 | 16 combinations 
5 | 25 combinations 
6 | 36 combinations 
Total | 91 combinations 


The relationship between the number of combinations is determined by taking a sum of the square value of the dice. 


6 
; n® 
n=1 


Where n is the value the dice can be, and A» is the number of dice Team B rolls. 


This needs to be divided by the number of total possible combinations: 
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Team A | Team B 

1 | (1 and 1) or (1 and 2) or (1 and 3) etc.. 
2 | (1 and 1) or (1 and 2) or (1 and 3) etc.. 
etc... 


The relationship to determine the number of total combinations is determined by taking 6 to the power of the number of dice from boths teams: 
6417 A2 


Therefore, the probability of Team A winning if it only rolls one dice is: 


6 
aaa 
Ai +A 
6 1 2 nel 


However, we need to account for Team A having more than one dice. Imagine Team A has two dice, and Team B has one dice, which times does 


Team A win: 
Team A Team B 
1 and other (1) and (1) 
2 and other (1 or 2) and (1 or 2) 
or 2 or 3) and (1 or 2 or 3) 


(1 or 2 or 3 or 4) and (1 or 2 or 3 or 4) 


w 

@ 
3 
a 
fo} 
fa 
> 
oO 
3 

=~ 
b 


2 
(1 or 2 or 3 or 4 or 5) and (1 or 2 or 3 or 4 or 5) 
(1 or 2 or 3 or 4 or 5 or 6) and (1 or 2 or 3 or 4 or 5 or 6) 
Becuase Team A wins if the first dice is the highest, the second dice is either equal to or less than the first dice, therefore, the actual amount of 
times that Team A wins is evry combination of these two dice: 
Team A | Team B 
1 and 1 | (1) and (1) 
(2 and 1) or (2 and 2) or (1 and 2) | (1 or 2) and (1 or 2) 


etc. 


Count the number of different combinations that can occur: 


Team A | Team B 
1 combination | (1) and (1) 
3 combinations | (1 or 2) and (1 or 2) 
5 combinations | (1 or 2 or 3) and (1 or 2 or 3) 
7 combinations | (1 or 2 or 3 or 4) and (1 or 2 or 3 or 4) 
9 combinations | 
| 


2 
(1 or 2 or 3 or 4 or 5) and (1 or 2 or 3 or 4 or 5) 
2 


11 combinations (1 or 2 or 3 or 4 or 5 or 6) and (1 or 2 or 3 or 4 or 5 or 6) 


The number of combinations therefore is: 


ni — (n — 1): 


Multiply the times Team A can win by all the combinations Team A can win, you get the Probability of Team A winning Phase 1: 


1 6 
P() Phase = gata ) Ae —(n- 1)*") an” 
n= 


and the probability of Team A not winning Phase 1: 


lis P() Phaset 


#Probability of Team A winning Phase 1 
def p_phase1(team_A_count, team_B_count): 
#Different values of the dice 
ns = [1, 2, 3, 4, 5, 6] 


#Calculate the sum of (n*{A_1}-(n-1)*{A_1}) * n*{A_2} between 1 and 6 
_sum = sum(((n**team_A_count) - ((n - 1)**team_A_count)) * (n**team_B_count) for n in ns) 


#Calculate the total number of combinations \frac{1}{6*{A_1+A_2}} 
quantity = 1/(6**(team_A_count+team_B_count)) 
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#Return the product of these two parts \frac{1}{6*{A_1+A_2}} \sum_{n=1}*{6} (n*{A_1}-(n-1)*{A_1}) * n*{A_2} 
return quantity * _sum 


Item 2: Determine the probability of Team A wins Phase 1 in k amount of rounds out of j total 


Vv 


rounds. 


¢ Out of j amount of turns, how many times does Team A win Phase 1. 
Consider, S=Success and F=Failure, where Success=Team A Winning and Failure is Team A not Winning. 


Out of 10 turns, the following could occur: 


Turn 112345678910 
Outcome |S FFFFSFFFS 


This is 3 Successes, and therefore 7 Failures. 
Probability of Winning 3 Phase 1's in 10 Turns then is: 
P()i/j =ps*pr*pr* pr * pr * ps * pr * pr * Pr * ps 
Where c=Number of times Team A win Phase 1, and j is the total number of turns conducted, Pg is probability of success, and Pr is probability 
of failure. Where k=3 (number of success), and j=10 (total number). 


This can be simplified to: 


PQ k/j = ps® * pp?) 


This can actually be made agnostic to the actual numbers: 


P()k/j = ps* « ppo*) 


So far though, this is just the probability of that exact scenario above, there is actually multiple ways of getting 3 wins out of 10: 


Outcome Variation 1 | SFFFFSFFFS 
Outcome Variation 2 | FSFFFSFFFS 


The number of combinations that are available is determined by a factorials: 
j! 
k! x I! 
Where | is the amount of failures. 


The probability then is: 
! k I 
PTeamAwinningkofj = Kall *DPs * PF 


Where pr = 1 — pg, therfore: 


PTeamAwinningkofj = klall * ps* * (1 — ps)! 


and further simplfications: 
i! 


PTeamAwinningkofj — klix (j—h)! * ps* * (1 - ps)0*) 


This is known as a binomial distribution. 


We can see what the probability of winning j turns out of 10 total turns, by setting j=10, and k=1, 2, 3, 4 etc. 


import math 
import numpy as np 
import matplotlib.pyplot as plt 


def p_TeamA_Winning_k_of_j_turns(team_A_count, team_B_count, win_count, total_count): 
#Calculate the quantity of combinations 
j = total_count 
k = win_count 
quantity = math.factorial(j) / ((math.factorial(k)) * (math. factorial(j-k))) 


#calculate the quantity of events of success and failure 
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prob = p_phase1(team_A_count, team_B_count)**k * (1-p_phase1(team_A_count, team_B count) )**(j-k) 


# Return probability 
return quantity * prob 


#Define count of dice on each side 
team_A_count = 4 
team_B_count = 8 


#Define the number of turns that will occur 
total_turn_count = 10 


#Go through each turn, from 1 to total turn count and capture the probability at each stage 
probs = [] 
for t in range(1, total_turn_count+1): 
probability = p_TeamA_Winning_k_of_j_turns(team_A_count, team_B_count, t, total_turn_count) 
probs. append(probability) 
#print(f"{int(t)} Round(s) out of {total_turn_count}: {probability}") 


#Plot each of the probabilities 

plt.figure(figsize=(10, 6)) 

bars = plt.bar(range(1, total_turn_count+1), probs, color='skyblue' ) 
plt.xlabel('Number of Wins for Team A in 10 Rounds - mathematical’) 
plt.ylabel( ‘Probability’ ) 

plt.title("Figure 1"); 
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This shows the probability of winning j rounds out of 10. 


¢ 1 Round(s) out of 10: 0.0019 
¢ 2 Round(s) out of 10: 0.0126 
¢ 3Round(s) out of 10: 0.0483 
¢ 4Round(s) out of 10: 0.1214 
e 5 Round(s) out of 10: 0.2091 
¢ 6 Round(s) out of 10: 0.2502 
¢ 7 Round(s) out of 10: 0.2052 
¢ 8 Round(s) out of 10: 0.1105 
¢ 9 Round(s) out of 10: 0.0352 
¢ 10 Round(s) out of 10: 0.0050 


The sum of these values is 1, which acts as a check to see if everything is working correctly. 


The second check we can do is simulating this problem, the following code simulates Team A and Team B, and keeps track of the amount of 
times Team A wins. 
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import random 


def simulate_Phase1_win_rate(team_A_count, team_B_count, total_count): 
# Initialise the number of wins of Team A 
team_A_Phasel_win = @ 


#Go through every turn 
for _ in range(total_count): 
# Team A rolls 


roll_a = max(random.randint(1, 6) for _ in range(team_A_count)) 


# Team B rolls 
roll_b = max(random.randint(1, 6) for _ in range(team_B_count) ) 


# Team A wins or it's a draw 
if roll_a >= roll_b: 
team_A_Phasel_win += 1 


#Return number of times team a won 
return team_A_Phase1_win 


#Game parameters 
total_turns = 10 
team_A_count = 4 
team_B_count = 8 


# Simulation parameters 
num_simulations = 10000 


# Run the Monte Carlo simulation 
team_A_Phase1_wins = [simulate_Phase1_win_rate(team_A_count, team_B_count, total_turns) for _ in range(num_simulations) ] 


# Plotting the probabilities as a histogram 

plt.figure(figsize=(10, 6)) 

plt.hist(team_A Phasel_wins, color='skyblue', alpha=0.7, density=True) 
plt.xlabel( ‘Number of Team A Wins in 10 Rounds - simulation’) 
plt.ylabel('Probability' ) 

plt.title("Figure 2"); 
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This shows the probability of winning j rounds out of 10. 


* 1 Round(s) / 10: 0.0019 
* 2Round(s) / 10: 0.0115 
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¢ 3 Round(s) / 10: 0.048 

¢ 4Round(s) / 10: 0.1233 
e 5 Round(s) / 10: 0.2079 
¢ 6 Round(s) / 10: 0.2443 
e 7 Round(s) / 10: 0.2099 
* 8 Round(s) / 10: 0.1136 
9 Round(s) / 10: 0.0345 
10 Round(s) / 10: 0.005 


This shows a very similar result to the results shown in Figure 1: 


¢ Simulation | Mathematical 

¢ 1 Round(s) / 10: 0.0019 | 0.0019 
¢ 2Round(s) / 10: 0.0115 | 0.0126 
¢ 3Round(s) / 10: 0.0480 | 0.0483 
¢ 4Round(s) / 10: 0.1233 | 0.1214 
e 5 Round(s) / 10: 0.2079 | 0.2091 
¢ 6 Round(s) / 10: 0.2443 | 0.2502 
¢ 7 Round(s) / 10: 0.2099 | 0.2052 
¢ 8 Round(s) / 10: 0.1136 | 0.1105 
¢ 9 Round(s) / 10: 0.0345 | 0.0352 
¢ 10 Round(s) / 10: 0.005 | 0.0050 


v |tem 3: Determine the probability of Taking 4 Health off Team A with no roll save. 


The problem with determining the probability of taking health off of Team A is the difference in probability between when not able to use roll 
saves, and using roll saves. For now, i consider these as two stages, no-roll-save and roll-save. 


Item 3 looks at no-roll-save 
Team A has 5 health and can only do a roll save when its health is at 1, therefore, what is the probabiltiy of losing 4 health in k amount of turns. 


To remove health from Team A, Team B rolls 16 (8x2) amount of dice, if any of the dice is a 6, then that dice is re-rolled, if it is 6 again then 
health is removed. Therefore the probability of success is: 


Therefore, probability of failure is: 
35 
36 


So what is the probability of f amount of units suceed in their attacks. It can again be explained with an binomial distribution as used in 


pl) neusand = (1 — Dl) wena = 


p_{TeamA winning k of j}: 


! 
PT B fulattacks fof. se gel (pe) 
‘ecamBsucces fulattacks fofg 7! g = A! 


def binomial_probability(p_success, p_failure, total_attemps, total_successes): 
#define parameters 
g = total_attemps 
f = total_successes 


#binomial 

binom = math. factorial(g)/(math.factorial(f)*math. factorial (g-f)) 
succ = p_success**f 

fail = p_failure**(g-f) 


#Total probabiltiy 
return binom * succ * fail 


There are many different ways of wounding a unit 4 times in r amount of turns, for example: 


1FFRF2FEF4 
FRFF1111F 
FRFF4FRFEF 
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Fa 


ure (no succesfull attacks), and 1 and 2 are amout of succesfulll attacks. 


This is a variation of binomial probability, calls the multinomial probability. The equatio for a multinomial probability is: 


n! 
P(X, = @1,.....,X_ = 2k) = 7 * (pi)... (pr) 
L1+..- LE: 
So the probability of getting 1 or 2 or 3 or 4 wounds in r turns would be: 
rl 
P(X), Xo, X3,.X4) = sae * Pl & Py? & PZ? « Pot x Py 


Where: 


X_1 = Number 
X_2 = Number 
X_3 = Number 
X_4 = Number 
X_F = Number 


P_1 = p(Team 
P_2 = p(Team 
P_3 = p(Team 
P_4 = p(Team 
1 - (PL 


ao) 
1 
mn 
" 


Now, determine 


succesfull 
succesfull 
succesfull 


succesfull 


PF NN WwW SB 


succesfull 


= Lei kalke! 


of 1 succesfull attacks 
of 2 succesfull attacks 
of 3 succesfull attacks 
of 4 succesfull attacks 


of F succesfull attacks 


winning Phase 1 j out of k times, and 1 succefull attack) = p_{phase1TeamBkofj}*p_{TeamB succesful attacks 1 of g} 
winning Phase 1 j out of k times, and 2 succefull attacks) = p_{phase1TeamBkofj}*p_{TeamB succesful attacks 2 of g} 


B 
B 
B winning Phase 1 j out of k times, and 3 succefull attacks) = p_{phaselTeamBkofj}*p_{TeamB succesful attacks 3 of g} 
B winning Phase 1 j out of k times, and 4 succefull attacks) = p_{phaselTeamBkofj}*p_{TeamB succesful attacks 4 of g} 
1 


+ P_2 + P_3 + P_4) 
all the combinations of 1, 2, 3, 4 which make 4 sccessfull attacks out of r turns: 


attacks in one round 

attacks in one round, 1 succesfull attack in one round 

attacks in one round, 2 succesfull attacks in one round 

attacks in one round, 1 succesfull attack in one round, 1 succesfull attack in one round 


attack in one round, 1 succesfull attack in one round, 1 succesfull attack in one round, 1 succesfull attack in one round 


To determine the total probability, sum P(X_1, X_2, X_3, X_4) for each combination: 


def p_TeamB_Win 
#Calculate th 
j = total_cou 
k = win_count 
quantity = ma 


#calculate th 


rl 
oy * P®l x P® x P®™ ~ P™ % Pet 
1 2 3 4 F 
ae{(4),(3+1),(2+2),(2+1-+1),(1+1+141)} Xy!Xy1X3! Kal Xr! 


ning _k_of_j_turns(team_A_count, team_B_ count, win_count, total_count): 
e quantity of combinations 

nt 

th.factorial(j) / ((math.factorial(k)) * (math. factorial(j-k))) 


e quantity of events of success and failure 


prob = (1-p_phase1(team_A_count, team_B_count))**k * p_phase1(team_A_count, team_B_count)**(j-k) 


# Return probability 


return quanti 


ty * prob 


def probability_of_TeamA_losing 4 health(r, team_A_count, team_B_count, p_attack_succesfull, p_attack_fail): 


#Write out al 


# Where; 

# [@, @, 
# [1, @, 
# [@, 2, 
# F255. als 
# [4, ®, 
combinations 


#Initialise t 
probability = 


1 the combinations 


@, 1] = 1 round of 4 succesfull attack (4 total successfull attacks) 

1, @] = 1 rounds of 1 succesfull attacks, 1 round of 3 succesfull attacks (4 total successfull attacks) 
@, @] = 2 rounds of 2 succesfull attacks (4 total successfull attacks) 

®, @] = 2 rounds of 1 succesfull attacks, 1 round of 2 succesfull attacks (4 total successfull attacks) 
®, @] = 4 rounds of 1 succesfull attacks (4 total successfull attacks) 


= [[@, @, @, 1], [1, @, 1, ®], [@, 2, @, @], [2, 1, @, 6], [4, @ @, @]] 


he probability 
) 


#Go through each combination and sum up the probability 


for combinati 
#Each combi 


on in combinations: 
nation has a different amount of succesfull wins i.e. 
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succesfull_attacks = sum(combination) 


if (succesfull_attacks > r): 
succesfull_attacks = r 


#Calculate the amount of failures. 
failed_attacks = r - succesfull_attacks 


#Calculate the probability of winning succesfull_attacks out of r amount of turns 
p_phaselTeamBkofj = p_TeamB_Winning_k_of_j_turns(team_A_count, team_B_count, succesfull_attacks, r) 


#Calculate the probability of getting 1, 2, 3, 4 succesfull attacks in 16 dice (Team B has 8 dice, multiply by 2 for 16) 
p_attacks1 = binomial_probability(p_attack_succesfull, p_attack_fail, team_B_count*2, 1) 
p_attacks2 = binomial_probability(p_attack_succesfull, p_attack_fail, team_B_count*2, 2) 
p_attacks3 = binomial_probability(p_attack_succesfull, p_attack_fail, team_B_count*2, 3) 
p_attacks4 = binomial_probability(p_attack_succesfull, p_attack_fail, team_B_count*2, 4) 


#Product of probabilities. 

p_1 = p_phaselTeamBkofj * p_attacks1 
p_2 = p_phaselTeamBkofj * p_attacks2 
p_3 = p_phaselTeamBkofj * p_attacks3 
p_4 = p_phaselTeamBkofj * p_attacks4 
pf=1-(pi+p2+p3+p 4) 


#Put everything together 

multinomial = math.factorial(r)/((math.factorial(combination[®])*math. factorial (combination[1])*math. factorial (combination[2])*math. fact 
prob = (p_1**combination[@]) * (p_2**combination[1]) * (p_3**combination[2]) * (p_4**combination[3]) * (p_f**failed_attacks) 

delta = (multinomial * prob) 


#Add the delta probability to the total sum of probabilities 
probability += delta 


#Return the probability 
return probability 


#Define count of dice on each side 
team_A_count = 4 
team_B_count = 8 


#Define the maximum number of turns to calculate probability over 
max_turns = 48 


#Determine the probability of one unit doing a succesfull attack (two d6 in a row) 
p_attack_succesfull = (1/6 * 1/6) 
p_attack_fail = 1 - p_attack_succesfull 


#Go through each turn up to the max and calculate the probability 

probs = [] 

for r in range(1, max_turns): 
prob = probability_of_TeamA_losing 4 health(r, team_A_count, team_B_count, p_attack_succesfull, p_attack_fail) 
probs. append(prob) 
#print(f"{t} Round(s) / {max_turns}: {prob}") 


#Sum the total probability (This should be equal to 1?) 
print (sum(probs) ) 


#Plot the histogram 

plt.figure(figsize=(10, 6)) 

bars = plt.bar(range(1, max_turns), probs, color='skyblue') 

plt.xlabel( ‘Probability of TeamA losing 4 health with no save rolls - mathematical’) 
plt.ylabel( ‘Probability’ ) 

plt.title("Figure 3"); 
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Run a simulation to test the probabilities against. 


def simulate_TeamB_Succesfully making 4 attacks(team_A_count, team_B count, save_roll=False): 
# Team A stats - health to remove 
health_to_remove_a = 4 


# Simulate rounds until Team A loses all health 


total_turns 


= @ 


#Continue until all 4 health is gone 

while health_to_remove_a > @: 
#Increment turn count 
total_turns += 1 


# Team A 
roll_a = 


# Team B 
roll_b = 


rolls 4 dice 


max(random.randint(1, 6) for _ in range(team_A_count)) 


rolls 8 dice 


max(random.randint(1, 6) for _ in range(team_B_count)) 


if roll_a >= roll_b: # Team A wins or it's a draw 
continue 


else: # 


Team B wins 


#Determine health loss 


for _ 


in range(team_B_count*2): 


die = random.randint(1, 6) #Roll dice 


return total_ 


#Define count of 
team_A_count = 4 
team_B_count 8 


if die == 6: # If dice shows a 6, re-roll 


re_roll = random.randint(1, 6) 
if re_roll == 6: # If the re-roll is also a 6, Team A loses health 
if (save_roll): 
_save_roll = random.randint(1, 6) 
if (_save_roll == 1): 
health_to_remove_a -= 1 
else: 
health_to_remove_a -= 1 
turns 


dice on each side 


# Simulation parameters 
num_simulations = 10000 
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# Run the Monte Carlo simulation 
turn_counts = [simulate_TeamB_Succesfully_making_4 attacks(team_A_count, team_B_count) for _ in range(num_simulations) ] 


# Plotting the probabilities as a histogram 

plt.figure(figsize=(10, 6)) 

plt.hist(turn_counts, bins=20@, color='skyblue', alpha=@.7, density=True) 
plt.xlabel('Probability of TeamA losing 4 health with no save rolls - simulation') 
plt.ylabel( ‘Probability’ ) 

plt.title("Figure 4"); 


(4 


Figure 4 


0.07 


0.06 


0.05 


0.04 


Probability 


0.03 


0.02 


0.01 


0.00 T 
0 20 40 60 80 100 


Drohahility af TaamA Incinn 4 health with na cave roiic - cimilatinn 
This is the time when things stop matching. 


v Item 4: Determine the probability of Taking 4 Health off Team A with roll saves. 


Now work out the probability when save rolls are allowed. 


A save roll works, when Team A has 1 health left, and Team B makes a succesfull attack (2 d6 in a row) then Team B can roll one dice, on a 2 or 


above the succesfull attack is void. However, Team A has 3 skill points that it can use to upgrade a 1 to a 2. 
Therefore, essentially, after the initial 4 health is removed from Team A, then Team B has to roll a 1 four times. 


So, the probability of Team B making a succesfull attack is now: 


1 1 1 
P()wound oe 


So technically, it should be the same analysis as before, except now the probability of succesfull attacks is higher 


#Define count of dice on each side 
team_A_count = 4 
team_B_count = 8 


#Define the maximum number of turns to calculate probability over 
max_turns = 600 


#Determine the probability of one unit doing a succesfull attack (two d6 and 1 di in a row) 
p_attack_succesfull = (1/6 * 1/6 * 1/6) 
p_attack_fail = 1 - p_attack_succesfull 


#Go through each turn up to the max and calculate the probability 
probs = [] 
for r in range(1, max_turns): 
prob = probability_of_TeamA_losing_4 health(r, team_A_count, team_B_count, p_attack_succesfull, p_attack_fail) 


probs. append(prob) 
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#print(f"{t} Round(s) / {max_turns}: {prob}") 


#Sum the total probability (This should be equal to 1?) 
print (sum(probs) ) 


#Plot the histogram 

plt.figure(figsize=(10, 6)) 

bars = plt.bar(range(1, max_turns), probs, color='skyblue' ) 

plt.xlabel( ‘Probability of TeamA losing 4 health with save rolls - mathematical’) 
plt.ylabel( ‘Probability’ ) 

plt.title("Figure 5"); 
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This isnt right atall, value is further away from 1.0. 


Compare it to a simulation. 


#Define count of dice on each side 
team_A_count = 4 
team_B_count = 8 


# Simulation parameters 
num_simulations = 5000 


# Run the Monte Carlo simulation 
turn_counts = [simulate_TeamB_Succesfully_making_4 attacks(team_A_count, team_B_count, save_roll=True) for _ in range(num_simulations) ] 


#Print all probabilities 
#for i, t in enumerate(np.arange(1, total_turns + 1)): 
#print(f"{t} Round(s) / 10: {probabilities[i]}") 


# Plotting the probabilities as a histogram 

plt.figure(figsize=(10, 6)) 

plt.hist(turn_counts, bins=20@, color='skyblue’, alpha=0@.7, density=True) 
plt.xlabel( ‘Probability of TeamA losing 4 health with no save rolls - simulation’) 
plt.ylabel( ‘Probability’ ) 

plt.title("Figure 6"); 
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Mathematical does not match simulation at all. 


y_ Item 5: Determine the probability of taking all Health off off Team A. 


Now Probabilities of taking 4 health off with no save rolls, and taking health off with save rolls, i should just be able to add these together 
somehow? Is that correct? It didnt work. 
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#Define count of dice on each side 
team_A_count = 4 
team_B_count = 8 


#Define the maximum number of turns to calculate probability over 
max_turns = 180 


#Determine the probability of one unit doing a succesfull attack (two d6 in a row) 
p_attack_succesfull = (1/6 * 1/6) 
p_attack_fail = 1 - p_attack_succesfull 


#Determine the probability of one unit doing a succesfull attack (two d6 and 1 d1 in a row) 
p_attack_succesfull_with_save_roll = (1/6 * 1/6 * 1/6) 
p_attack_fail_with_save_roll = 1 - p_attack_succesfull_with_save_roll 


#Go through each turn up to the max and calculate the probability 
probs = [] 
for r in range(1, max_turns): 
#Calculate probability of with no save roll 
prob = probability_of_TeamA_losing 4 health(r, team_A_count, team_B_count, p_attack_succesfull, p_attack_fail) 


#Add probability of with save roll 
prob += probability_of_TeamA_losing 4 health(r, team_A_count, team_B_count, p_attack_succesfull_with_save_roll, p_attack_fail_with_save_rc 


#Add to list to plot 
probs. append (prob) 


#Sum the total probability (This should be equal to 1?) 

print(sum(probs) ) 

#Plot the histogram 

plt.figure(figsize=(10, 6)) 

bars = plt.bar(range(1, max_turns), probs, color="skyblue') 
plt.xlabel('Probability of TeamA losing 4 health with save rolls - mathematical') 
plt.ylabel( ‘Probability’ ) 

plt.title("Figure 7"); 
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This is definitly also incorrect. 


And | did a simulation of the whole simulation to make sure. We also get the probability of the numebr of Team B deaths aswell. 


def Full_Game(team_A_count, team_B_count): 
# Team A stats - health to remove 
health_to_remove_a = 5 
skill_points_a = 3 
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# Simulate rounds until Team A loses all health 
total_turns = 1 

teamB_killed = @ 

succesfull _save_rolls = @ 


#Continue until all 4 health is gone 
while health_to_remove_a > @: 
# Team A rolls 4 dice 


roll_a = max(random.randint(1, 6) for _ in range(team_A_count) ) 


# Team B rolls 8 dice 


roll_b = max(random.randint(1, 6) for _ in range(team_B_count)) 


if roll_a >= roll_b: # Team A wins or it's a draw 


#roll team_B_count dice, if any dice is over 2 then a unit of team B is killed 


for _ in range(team_B_count): 
dice = random.randint(1, 6) 
if dice > 2: 
teamB_killed += 1 
else: # Team B wins 
#Determine health loss 
for _ in range(team_B_count*2): 
die = random.randint(1, 6) #Roll dice 
if die == 6: # If dice shows a 6, re-roll 
re_roll = random.randint(1, 6) 
if re_roll == 6: # If the re-roll is also a 6, Team A loses health 
#If potential health loss would kill Team A then do save roll 
if health_to_remove_a - 1 > @: 
health_to_remove_a -= 1 
else: 
#else do a save roll 
save_roll = random.randint(1, 6) 


#If save roll greater than 1 then dont remove health, if it is 1 then use skill point to save (if points left over) 


if save_roll == 
if skill_points_a == 
health_to_remove_a -= 1 
else: 
#Use skill_points_a to save 
skill_points_a -= 1 
succesfull_save_rolls += 1 
else: 
succesfull_save_rolls += 1 


#Increment turn count 
total_turns += 1 


return total_turns, teamB_killed, succesfull_save_rolls 


#Define count of dice on each side 
team_A_count = 4 
team_B_count = 8 


# Simulation parameters 
num_simulations = 50e0 


# Run the Monte Carlo simulation 

# Initialize empty lists to store results 
total_turns = [] 

teamB_killeds = [] 

succesfull_save_rolls = [] 


# Run the simulations and unpack each result 
for _ in range(num_simulations): 
total_turn, teamB_killed, save_rolls = Full_Game(team_A_count, team_B_count) 
total_turns.append(total_turn) 
teamB_killeds.append(teamB_killed) 
succesfull_save_rolls.append(save_rolls) 


# Set up the figure and axes 
plt.figure(figsize=(15, 5)) 


#Probability of Reaching round n 

plt.subplot(1, 3, 1) 

plt.hist(total_turns, bins=200, color='skyblue', alpha=@.7, density=True) 
plt.title("Number of Rounds Reached") 

plt.xlabel("Turn Count") 
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plt.ylabel("Probability") 


#Probability of Team B Deaths 

plt.subplot(1, 3, 2) 

plt.hist(teamB_killeds, bins=20@, color='skyblue', alpha=@.7, density=True) 
plt.title("Number of Team B Killed") 

plt.xlabel("Turn Count") 

plt.ylabel("Probability") 


#Probability of Save Rolls 

plt.subplot(1, 3, 3) 

plt.hist(succesfull_save_rolls, bins=200, color='skyblue', alpha=@.7, density=True) 
plt.title("Number of Save Rolls Made") 

plt.xlabel("Turn Count") 

plt.ylabel("Probability") 


# Show the plots 
plt.tight_layout() # Adjust layout to prevent overlap 


plt.show() 
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v_ Item 6: Determine the probability of killing j amount of Team B units. 


Wont do this yet, but i think this is easy, just a binomial thing? 


def p_TeamA_Killing_k_of_j_TeamBUnits(success_count, total_count): 
#Calculate the quantity of combinations 
j = total_count 
k = success_count 
quantity = math.factorial(j) / ((math.factorial(k)) * (math.factorial(j-k))) 


#Binomial of sucesses 
probability_success = binomial_probability(4/6, (1-4/6), 8, success_count) 


#calculate the quantity of events of success and failure 
prob = probability_success**k * (1-probability_success)**(j-k) 


# Return probability 
return quantity * prob 


#Define the number of turns that will occur 
teamA_Count = 8 


#Go through each turn, from 1 to total turn count and capture the probability at each stag 
probs = [] 
for t in range(1, teamA_Count): 

probability = binomial_probability(4/6, (1-(4/6)), teamA_Count, t) 

probs. append(probability) 
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#print(f"{int(t)} Round(s) out of {total_turn_count}: {probability}") 


#Plot each of the probabilities 

plt.figure(figsize=(10, 6)) 

bars = plt.bar(range(1, teamA_Count), probs, color='skyblue' ) 
plt.xlabel('Number of Wins for Team A in 10 Rounds - mathematical’) 
plt.ylabel('Probability' ) 

plt.title("Figure 1"); 


sum(probs ) 
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Probably amount of rounds won out of Conquest Creatiosn 157 


#Define count of dice on each side 
team_A_count = 4 
team_B_count = 8 


#Define the number of turns that will occur 
total_turn_count = 157 


#Go through each turn, from 1 to total turn count and capture the probability at each stage 
probs = [] 
for t in range(1, total_turn_count+1): 
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